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ABSTRACT 

I review recent progress in understanding the formation of the first stars and quasars. The 
initial conditions for their emergence are given by the now firmly established model of cosmological 
structure formation. Numerical simulations of the collapse and fragmentation of primordial gas 
indicate that the first stars formed at redshifts z ~ 20 — 30, and that they were predominantly 
very massive, with M* > 100M Q . Important uncertainties, however, remain. Paramount among 
them is the accretion process, which builds up the final stellar mass by incorporating part of the 
diffuse, dust-free envelope into the central protostellar core. The first quasars, on the other hand, 
are predicted to have formed later on, at z ~ 10, in more massive dark matter halos, with total 
masses, ~ 10 8 M©, characteristic of dwarf galaxies. 

Subject headings: cosmology: theory — galaxies: high-redshift — stars: formation 



1. INTRODUCTION 

How and when did the first sources of light 
form in the universe? Within the framework of 
modern cosmology, we have learned that the first 
stars and quasars formed at the end of the so- 
called "dark ages" , a few 10 s yr after the big 
bang (e.g., Barkana & Loeb 2001; Miralda-Escude 
2003). The cosmic dark ages began ~ 500,000 yr 
after the big bang when the photons of the cos- 
mic microwave background (CMB) were emitted. 
At this time, the CMB photons redshifted into 
the infrared wavelength band so that they were no 
longer energetic enough to ionize hydrogen atoms, 
and could henceforth propagate freely. From the 
viewpoint of a human observer, the universe de- 
scended into a state of complete darkness, that 
lasted until the cosmic Renaissance of first light 
occurred much later on. 

At these early times, the universe did not yet 
contain any gravitationally bound structures, and 
its state was simple enough to allow a descrip- 
tion in terms of precise, linear physics (e.g., Hu 
& Dodelson 2002). We know from observations 
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of the highest redshift quasars and galaxies that 
the dark ages must have ended prior to z ~ 7 
(e.g., Hu et al. 2002; Fan et al. 2003). Whereas 
we can directly probe the state of the universe at 
z ~ 1000, by measuring CMB anisotropies, and at 
z ^5 6, by observing quasar absorption lines and 
Lyman-a emission from galaxies, we know little 
about the formative epoch in between. Thus, this 
era is the final frontier of observational cosmol- 
ogy, and the crucial question is (e.g., Barkana & 
Loeb 2001; Loeb & Barkana 2001; Miralda-Escude 
2003): How and when did the dark ages end? 

Within the now firmly established model of cos- 
mological structure formation (e.g., Spergel et al. 
2003) , the first stars are predicted to have formed 
at redshifts z ~ 20 — 30 in dark matter (DM) ha- 
los of total mass ~ 10 6 Af Q (e.g., Couchman & 
Rees 1986; Ostriker & Gnedin 1996; Tegmark et al. 
1997; Abel et al. 1998; Yoshida et al. 2003a). Nu- 
merical simulations of the collapse and fragmen- 
tation of metal-free gas, assumed to contain no 
dynamically significant magnetic fields yet, have 
indicated that the first stars, the so-called Popu- 
lation III (e.g., Bond 1981), were predominantly 
very massive, with > IOOMq (e.g., Bromm, 
Coppi, & Larson 1999, 2002; Abel, Bryan, & Nor- 
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man 2000, 2002; Nakamura & Umemura 2001). 
Many uncertainties, however, remain. Two key 
problems are: (i) The physics of accretion from a 
dust-free envelope, which determines the final stel- 
lar mass, is not yet well understood (e.g., Omukai 
& Palla 2001, 2003; Tan & McKee 2003). (ii) 
A complete theory of primordial star formation 
would allow us to predict the exact functional 
form of the Population III initial mass function 
(IMF). Again, a reliable determination is still elud- 
ing us (see, e.g., Nakamura & Umemura 2001, 
2002; Omukai & Yoshii 2003 for first attempts). 
I will address both issues in §3. 

The first quasars, on the other hand, pow- 
ered by accretion onto supermassive black holes 
(SMBHs), are expected to have formed later on, at 
z ~ 10, in more massive, dwarf-sized systems (e.g., 
Haiman & Loeb 2001; Bromm & Loeb 2003a). Nu- 
merical simulations of high-redshift SMBH forma- 
tion indicate that these could have formed only 
after a significant episode of previous star forma- 
tion (e.g., Bromm & Loeb 2003a; Di Matteo et 
al. 2003). The long-standing question (e.g., Silk 
& Rees 1998) of what came first, stars or quasars, 
would then be answered in favor of the former. 

The formation of the first luminous sources had 
important implications for the evolution of the in- 
tergalactic medium (IGM). Massive Population III 
stars were efficient producers of ionizing photons 
(e.g., Tumlinson & Shull 2000; Bromm, Kudritzki, 
& Loeb 2001b; Schaerer 2002, 2003; Tumlinson, 
Shull, & Vcnkatesan 2003). The contribution from 
these stars, thus, could have been important in 
reionizing the IGM early on, at z > 15, as may 
be required by the recent Wilkinson Microwave 
Anisotropy Probe (WMAP) data (Kogut et al. 
2003; Spergel et al. 2003) on the optical depth 
to Thomson scattering (e.g., Cen 2003a,b; Cia- 
rdi, Fcrrara, & White 2003; Haiman & Holder 
2003; Holder et al. 2003; Kaplinghat et al. 2003; 
Sokasian et al. 2003a,b; Wyithe & Loeb 2003a,b). 
The second key feedback effect concerns the pro- 
duction of the first heavy elements, and their sub- 
sequent dispersal by supernova (SN) explosions 
into the pristine IGM (e.g., Gnedin & Ostriker 
1997; Madau, Ferrara, & Rees 2001; Mori, Ferrara, 
& Madau 2002; Bromm, Yoshida, & Hernquist 
2003; Wada & Venkatesan 2003). Once the IGM 
has been metal-enriched above a minimum level, 
star formation is predicted to shift from a predom- 



inantly high-mass (Population III) to a lower- mass 
(Population II) mode (e.g., Omukai 2000; Bromm 
et al. 2001a; Bromm & Loeb 2003b; Schneider et 
al. 2003a). The nature of this fundamental tran- 
sition, whether gradual or sudden, will depend on 
the detailed IGM enrichment history (e.g., Schnei- 
der et al. 2002; Mackey, Bromm, & Hernquist 
2003; Scannapieco, Schneider, & Ferrara 2003; Ri- 
cotti & Ostriker 2004; Yoshida, Bromm, & Hern- 
quist 2004). 

Very recently, novel empirical probes have be- 
gun to test our theoretical ideas on star and 
galaxy formation in the early universe. Next 
to the WMAP data on the ionization state of 
the IGM, these probes include determinations 
of the chemical abundance pattern in extremely 
metal-poor Galactic halo stars (e.g., Christlieb 
et al. 2002), measurements of the near IR cos- 
mic background radiation (e.g., Santos, Bromm, 
& Kamionkowski 2002; Salvaterra & Ferrara 
2003), and the prospect of detecting high-redshift 
gamma-ray bursts (GRBs) with the Swift satellite 
to be launched in 2004 (e.g., Lamb & Reichart 
2000; Bromm & Loeb 2002). 

In this review, I will briefly summarize the cur- 
rent state of key aspects in this rapidly evolving 
field (see Bromm & Larson 2004 for a more de- 
tailed discussion). 

2. BASIC COSMOLOGICAL FRAME- 
WORK 

The establishment of the current standard 
ACDM model for cosmological structure forma- 
tion has provided a firm framework for the study 
of the first stars (Spergel et al. 2003). Within 
variants of the CDM model, where larger struc- 
tures are assembled hierarchically through succes- 
sive mergers of smaller building-blocks, the first 
stars are predicted to form in DM minihalos of 
typical mass ~ 1O 6 M at redshifts z ~ 20 — 30 
(e.g., Couchman & Rees 1986). The virial temper- 
atures in these low- mass halos, T v ; r oc M 2 / 3 (l + z) 
(Barkana & Loeb 2001), are below the threshold, 
~ 10 4 K, for efficient cooling due to atomic hy- 
drogen (e.g., Oh & Haiman 2002). It was realized 
early on that cooling in the low-temperature pri- 
mordial gas had to rely on molecular hydrogen 
instead (Saslaw & Zipoy 1967; Peebles & Dickc 
1968). 
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Since the thermodynamic behavior of the pri- 
mordial gas thus is primarily controlled by H2 
cooling, it is crucial to understand the non- 
equilibrium chemistry of H2 formation and de- 
struction (e.g., Lepp & Shull 1984; Anninos & 
Norman 1996; Abel et al. 1997; Galli & Palla 
1998; Puy & Signore 1999). In the absence of 
dust grains to facilitate their formation (e.g., Hi- 
rashita & Ferrara 2002), molecules have to form 
in the gas phase. The most important formation 
channel is given by the sequence: H + e~ — > 
H~ +7, followed by H~ + H -> H 2 + e~ (McDow- 
ell 1961). The free electrons act as catalysts, and 
are present as residue from the epoch of recombi- 
nation (Seager, Sasselov, & Scott 2000), or result 
from collisional ionization in accretion shocks dur- 
ing the hierarchical build-up of galaxies (e.g., Mac 
Low & Shull 1986; Shapiro & Kang 1987). The 
formation of hydrogen molecules thus ceases when 
the free electrons have recombined. Calculations 
of H2 formation in collapsing top-hat overden- 
sities, idealizing the virialization of dark matter 
halos in CDM cosmogonies, have found a simple 
approximate relationship between the asymptotic 
H2 abundance and virial temperature in the over- 
density (or halo): /h 2 oc (Tegmark et al. 
1997). 

Applying the familiar criterion (Rees & Os- 
triker 1977; Silk 1977) for the formation of galaxies 
that the cooling timescale has to be shorter than 
the dynamical timescale, i coo i < tdyn, one can de- 
rive the minimum halo mass at a given redshift 
inside of which the gas is able to cool and eventu- 
ally form stars (e.g., Tegmark et al. 1997). The 
H2 cooling function has been quite uncertain over 
the relevant temperature regime, but recent ad- 
vances in the quantum-mechanical computation of 
the collisional excitation process (H atoms collid- 
ing with H2 molecules) have provided a much more 
reliable determination (see Galli & Palla 1998, and 
references therein). Combining the CDM prescrip- 
tion for the assembly of virialized DM halos with 
the H2 driven thermal evolution of the primordial 
gas, a minimum halo mass of ~ 10 6 M Q is required 
for collapse redshifts z vil . ~ 20 — 30. From de- 
tailed calculations, one finds that the gas in such 
a 'successful' halo has reached a molecule fraction 
in excess of /h 2 ~ 10~ 4 (e.g., Haiman, Thoul, & 
Loeb 1996; Tegmark et al. 1997; Yoshida et al. 
2003a). These systems correspond to 3 — 4<r peaks 




Fig. 1. — Emergence of primordial star forming re- 
gions within a standard ACDM cosmology. Shown 
is the projected gas density at z = 17 within a sim- 
ulation box of physical size ~ 50 kpc. The bright 
knots at the intersections of the filamentary net- 
work are the sites where the first stars form. (From 
Yoshida et al. 2003a.) 

in the Gaussian random field of primordial density 
fluctuations. In principle, DM halos that are suf- 
ficiently massive to harbor cold, dense gas clouds 
could form at higher redshifts, z v - lr > 40. Such 
systems, however, would correspond to extremely 
rare, high-er peaks in the Gaussian density field 
(e.g., Miralda-Escude 2003). 

To more realistically assess the formation of 
cold and dense star forming clouds in the earliest 
stages of cosmological structure formation, three- 
dimensional simulations of the combined evolution 
of the DM and gas are required within a cosmolog- 
ical set-up (e.g., Ostriker & Gnedin 1996; Gncdin 
& Ostriker 1997; Abel et al. 1998). These stud- 
ies confirmed the important role of H2 cooling in 
low-mass halos at high z. Recently, the problem of 
forming primordial gas clouds within a fully cos- 
mological context has been revisited with high nu- 
merical resolution (Yoshida et al. 2003a, b,c). The 
resulting gas density field is shown in Figure 1 for a 
standard ACDM cosmology at z = 17. The bright 
knots at the intersections of the filamentary net- 
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work are the star forming clouds, having individ- 
ual masses (DM and gas) of ~ 10 6 M Q . 

Whether a given DM halo successfully hosts a 
cold, dense (T < 0.5T vir ,n H > 5 x 10 2 cm" 3 ) gas 
cloud can nicely be understood with the Rees- 
Ostriker criterion, as can be seen in Figure 2. 
Yoshida et al. (2003a) derive a minimum collapse 
mass of M cr it — 7 x 10 5 M Q , with only a weak 
dependence on collapse redshift (see also Haiman 
ct al. 1996; Fuller & Couchman 2000; Machacek, 
Bryan, & Abel 2001). The dynamical heating ac- 
companying the merging of DM halos has an im- 
portant effect on the thermal and chemical evolu- 
tion of the gas. Clouds do not successfully cool if 
they experience too rapid a growth in mass. 

The primordial gas clouds that are found in the 
cosmological simulations are the sites where the 
first stars form. It is, therefore, important to learn 
what properties these clouds have in terms of over- 
all size, shape, and angular momentum content. 
The latter is often expressed by the familiar spin- 
parameter A = L\E\^ 2 /(GM 5 / 2 ), where L, E, 
and M are the total angular momentum, energy, 
and mass, respectively. The spin parameter is a 
measure of the degree of rotational support, such 
that the ratio of centrifugal to gravitational accel- 
eration is given by ~ A 2 at virialization. The spin 
values measured in pure DM cosmological simula- 
tions can be described by a lognormal distribution 
function with a mean value of A = 0.04, similar 
to what is found for larger-scale systems (Jang- 
Condell & Hernquist 2001). The overall sizes of 
the Population III star forming clouds are close 
to the virial radius of the host DM halo, with 
i? v i r ~ 100 pc, not too different from the typical 
dimensions of present-day giant molecular clouds 
(e.g., Larson 2003). Depending on the degree of 
spin, the clouds have shapes with various degree 
of flattening (e.g., Bromm et al. 2002; Yoshida 
et al. 2003a). Due to the importance of pressure 
forces, however, overall cloud shapes tend to be 
rather spherical. To fully elucidate the properties 
of the star-forming primordial clouds, even higher 
resolution cosmological simulations will be neces- 
sary. 

The theoretical predictions for the formation 
sites of the first stars sensitively depend on the 
exact nature of the DM component and its fluctu- 
ation spectrum. Recently, two models have been 
discussed that would much reduce the fluctua- 
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Fig. 2. — Condition for successful formation of 
cold, dense gas clouds. Mass-weighted mean H2 
fraction versus virial temperature. Open circles: 
DM halos that fail to host cold gas clouds. Filled 
circles: DM halos that succeed in harboring star- 
forming clouds. The critical line corresponds to 
the condition i cool ~ i ff . (From Yoshida et al. 
2003a.) 

tion power on small mass scales. The first of 
these, the warm dark matter (WDM) model (e.g., 
Bode, Ostriker, & Turok 2001), has been pro- 
posed to remedy the well-known problems of stan- 
dard ACDM on sub-galactic scales (e.g., Flores & 
Primack 1994; Moore et al. 1999). These con- 
cern the predicted large abundance of substructure 
or, equivalcntly, of satellite systems, and the high 
(cuspy) densities in the centers of galaxies. Both 
predictions are in conflict with observations. The 
second model, the 'running' spectral index (RSI) 
model, is suggested by the combined analysis of 
the WMAP data, the 2dF galaxy redshift survey, 
and Lyman-a forest observations (Spergel et al. 
2003; Peiris et al. 2003). A series of recent studies 
have worked out the consequences of these reduced 
small-scale power models on early star formation 
(Somerville, Bullock, & Livio 2003; Yoshida et al. 
2003b, c). Within these models, the star forma- 
tion rate at z > 15 is significantly reduced com- 
pared to the standard ACDM case. This is due to 
the absence of low-mass halos and their associated 
gas clouds that are cooled by molecular hydrogen. 
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Such a reduced rate of early star formation makes 
it difficult to achieve the large optical depth to 
Thomson scattering, as measured by WMAP, with 
reasonable choices for the star formation efficiency 
in the first galaxies. 

3. FORMING THE FIRST STARS 
3.1. Formation of Prestellar Clumps 

The metal-rich chemistry, magnetohydrody- 
namics, and radiative transfer involved in present- 
day star formation are complex, and we still lack 
a comprehensive theoretical framework that pre- 
dicts the IMF from first principles. Star formation 
in the high redshift universe, on the other hand, 
poses a theoretically more tractable problem due 
to a number of simplifying features, such as: (i) 
the initial absence of heavy elements and therefore 
of dust; (ii) the absence of dynamically-significant 
primordial magnetic fields; and (iii) the absence of 
any effects from previous episodes of star forma- 
tion which would completely alter the conditions 
for subsequent generations of stars. The cooling 
of the primordial gas depends only on hydrogen 
in its atomic and molecular form. Whereas in 
the present-day interstellar medium (ISM) the 
initial state of the star forming cloud is poorly 
constrained, the corresponding initial conditions 
for primordial star formation are simple, given by 
the popular ACDM model of cosmological struc- 
ture formation. 

How did the first stars form? This subject has 
a long and venerable history (e.g., Schwarzschild 
& Spitzer 1953; Yoncyama 1972; Hutchins 1976; 
Silk 1977, 1983; Yoshii & Sabano 1979; Carlberg 
1981; Kashlinsky & Rees 1983; Palla, Salpeter, & 
Stahler 1983; Yoshii & Saio 1986). In this review, 
I focus mainly on the more recent work since the 
renewed interest in high-rcdshift star formation 
beginning in the mid-1990s (e.g., Haiman ct al. 
1996; Uehara et al. 1996; Haiman & Loeb 1997; 
Tegmark et al. 1997; Larson 1998). The complete 
answer to this question would entail a theoreti- 
cal prediction for the Population III IMF, which 
is rather challenging. A more tractable task is to 
estimate the characteristic mass scale, M c , of the 
first stars, and most of the recent numerical work 
has focused on this simpler problem. The char- 
acteristic mass is the mass below which the IMF 
flattens or begins to decline (see Larson 1998 for 



examples of possible analytic forms). This mass 
scale is observed to be ~ 1M Q in the present-day 
universe. Since the detailed shape of the primor- 
dial IMF is highly uncertain, it is reasonable to 
first constrain M c , as this mass scale indicates the 
typical outcome of the primordial star formation 
process or, in other words, where most of the avail- 
able mass ends up (e.g., Clarke & Bromm 2003). 

To fully explore the dynamical, thermal, and 
chemical properties of primordial gas, three- 
dimensional numerical simulations are needed, al- 
though more idealized investigations in two, one, 
or even zero dimension are important in that they 
allow to probe a larger parameter space. The 
proper initial conditions for primordial star for- 
mation are given by the underlying model of cos- 
mological structure formation (see § 2). It is, 
therefore, necessary to simulate both the DM and 
gaseous ('baryonic') components. To date, two 
different numerical approaches have been used to 
simulate the general three-dimensional fragmenta- 
tion problem in its cosmological context. The first 
series of simulations used the grid-based adaptive 
mesh-refinement (AMR) technique (Abel ct al. 
2000, 2002; ABN henceforth), and the second one 
employed the smoothed particle hydrodynamics 
(SPH) method (Bromm et al. 1999, 2002; BCL 
henceforth). The SPH approach has the impor- 
tant advantage that it can easily accommodate the 
creation of sink particles (e.g., Bate, Bonnell, & 
Bromm 2003). Recently, the dynamical range of 
the standard SPH method has been significantly 
improved by implementing a 'particle splitting' 
technique (Kitsionas & Whitworth 2002; Bromm 
& Loeb 2003a). When the simulation reaches such 
high density in a certain region that a sink parti- 
cle would normally be created, a complementary 
strategy is now adopted: Every SPH particle in 
the unrefined, high-density region acts as a par- 
ent particle and spawns a given number of child 
particles, and endows them with its properties. In 
effect, this is the SPH equivalent of the grid-based 
AMR technique. 

The most important difference between these 
two studies lies in the way the initial conditions 
are implemented. The ABN simulations start at 
z = 100 with a periodic volume of physical size 
128/(1 +z) kpc. The AMR technique allows ABN 
to bridge the gap from cosmological to protostellar 
scales. The BCL effort, on the other hand, initial- 
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Fig. 3. — Properties of primordial gas. (a) Free 
electron abundance vs. hydrogen number density 
(in cm -3 ), (b) Hydrogen molecule abundance vs. 
number density, (c) Gas temperature vs. num- 
ber density. At densities below ~ 1 cm~ 3 , the gas 
temperature rises because of adiabatic compres- 
sion until it reaches the virial value of T vir ~ 5000 
K. At higher densities, cooling due to H2 drives 
the temperature down again, until the gas settles 
into a quasi- hydrostatic state at T ~ 200 K and 
n ~ 10 4 cm~ 3 . Upon further compression due 
to the onset of the gravitational instability, the 
temperature experiences a modest rise again, (d) 
Jeans mass (in M ) vs. number density. The 
Jeans mass reaches a value of Mj ~ 10 3 M Q for 
the quasi-hydrostatic gas in the center of the DM 
potential well. (From Bromm et al. 2002.) 

izes the simulations, also at z = 100, by realizing 
spherical overdensities that correspond to high-a 
peaks in the Gaussian random field of cosmological 
density fluctuations. 

In comparing the simulations of ABN and BCL, 
the most important aspect is that both studies, 
employing very different methods, agree on the ex- 
istence of a preferred state for the primordial gas, 
corresponding to characteristic values of tempera- 
ture and density, T c ~ 200 K, and n c ~ 10 4 cm~ 3 , 
respectively. These characteristic scales turn out 
to be rather robust in the sense that they are not 
very sensitive to variations in the initial condi- 
tions. In Figure 3 (panel c), this preferred state 
in the T — n phase diagram can clearly be dis- 



cerned. This figure, from the simulations of BCL, 
plots the respective gas properties for each individ- 
ual SPH particle. The diagram thus contains an 
additional dimension of information: where evolu- 
tionary timescales are short, only few particles are 
plotted, whereas they tend to accumulate where 
the overall evolution is slow. Such a 'loitering' 
state is reached at T c and n c . 

These characteristic scales can be understood 
by considering the microphysics of H2, the main 
coolant in metal-free, star forming gas (Abel et 
al. 2002; Bromm et al. 2002). At tempera- 
tures < 1000 K, cooling is due to the collisional 
excitation and subsequent radiative decay of ro- 
tational transitions. The two lowest-lying rota- 
tional energy levels in H2 have an energy spac- 
ing of Ej fee — 512K. Collisions with particles 
(mostly H atoms) that populate the high-energy 
tail of the Maxwell-Boltzmann velocity distribu- 
tion can lead to somewhat lower temperatures, 
but H 2 cooling cannot proceed to T < 100 K. 
This explains the characteristic temperature. The 
characteristic density, in turn, is given by the crit- 
ical density above which collisional de-excitations, 
which do not cool the gas, compete with radiative 
decays, which lead to cooling. This saturation of 
the H2 cooling marks the transition from NLTE ro- 
tational level populations to thermal (LTE) ones. 
At densities below n c , the cooling rate is propor- 
tional to the density squared, whereas at higher 
densities, the dependence is only linear. 

Once the characteristic state is reached, the 
evolution towards higher density is temporarily 
halted due to the now inefficient cooling, and the 
gas undergoes a phase of quasi-hydrostatic, slow 
contraction. To move away from this 'loitering' 
regime, enough mass has to accumulate to trigger 
a gravitational runaway collapse. This condition 
is simply M > Mj, where the Jeans mass or, al- 
most equivalently, the Bonnor-Ebert mass can be 
written as (e.g., Clarke & Bromm 2003) 

/ T \ 3/2 / n \-i/2 

Mj ~ 700M fl — ; 

J °V200 Ky V10 4 cm- 3 / 

(1) 

Primordial star formation, where magnetic fields 
and turbulence are expected to initially play no 
important dynamical role, may be the best case 
for the application of the classical Jeans criterion 
which is describing the balance between gravity 
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and the opposing thermal pressure. 

A prestellar clump of mass M > Mj is the im- 
mediate progenitor of a single star or, in case of 
further subfragmcntation, a binary or small multi- 
ple system. In Galactic star-forming regions, like 
p Ophiuchi, such clumps with masses close to stel- 
lar values have been observed as gravitationally 
bound clouds which lack the emission from em- 
bedded stellar sources (e.g., Motte, Andre, & Neri 
1998). The high-density clumps arc clearly not 
stars yet. To probe the further fate of a clump, 
one first has to follow the collapse to higher den- 
sities up to the formation of an optically thick hy- 
drostatic core in its center (see § 3.2), and sub- 
sequently the accretion from the diffuse envelope 
onto the central core (see § 3.3). The parent clump 
mass, however, already sets an upper limit for the 
final stellar mass whose precise value is determined 
by the accretion process. 

Although ABN and BCL agree on the mag- 
nitude of the characteristic fragmentation scale, 
ABN have argued that only one star forms per 
halo, while BCL have simulated cases where mul- 
tiple clumps form, such that the number of stars in 
a minihalo is A* ~ 1 — 5. Among the cases studied 
by BCL favoring N* > 1, are high-spin runs that 
lead to disk-like configurations which subsequently 
fragment into a number of clumps, or more mas- 
sive halos (with total mass M < 10 7 M Q ). Notice 
that the case simulated by ABN corresponds to a 
gas mass of ~ 3 x 10 4 M Q . Such a low-mass halo is 
only marginally able to cool, and when BCL simu- 
lated a system with this mass, they also found that 
only a single star forms inside the halo. Further 
work is required to elucidate whether the multiple- 
clump formation in BCL does also occur in simula- 
tions with realistic cosmological initial conditions. 

A possible clustered nature of forming Popula- 
tion III stars would have important consequences 
for the transport of angular momentum. In fact, a 
crucial question is how the primordial gas can so 
efficiently shed its angular momentum on the way 
to forming massive clumps. In the case of a clus- 
tered formation process, angular momentum could 
be efficiently transported outward by gravitational 
(tidal) torques. Tidal torques can then trans- 
fer much of the angular momentum from the gas 
around each forming clump to the orbital motion 
of the system, similar to the case of present-day 
star formation in a clustered environment (e.g., 



Larson 2002; Bate et al. 2003). In the isolated 
formation process of ABN, however, such a mech- 
anism is not available. Alternatively, ABN have 
suggested the efficient transport of angular mo- 
mentum via hydrodynamic shocks during turbu- 
lent collapse. Again, more work is required to con- 
vincingly sort out the angular momentum trans- 
port issue. 

Important contributions have also been made 
by more idealized, two- to zero-dimensional, stud- 
ies. These investigations typically ignore the dark 
matter component, thus implicitly assuming that 
the gas has already dissipatively collapsed to the 
point where the dark matter ceases to be dynam- 
ically dominant. Much attention has been paid 
to the collapse of filamentary clouds (Uehara et 
al. 1996; Nakamura & Umemura 1999, 2001, 
2002). Results from one- and two-dimensional 
simulations that otherwise include all the rele- 
vant processes for the primordial chemistry and 
cooling have estimated fragmentation scales M c ~ 
1M© and M c ~ 100M Q , depending on the (cen- 
tral) initial density. When sufficiently high (n > 
10 5 cm~ 3 ), the low-mass value is reached. This bi- 
furcation has led to the prediction of a bimodal 
IMF for the first stars (Nakamura & Umemura 
2001). It is however, not obvious how such a high 
initial density can be reached in a realistic situa- 
tion where the collapse starts from densities that 
are typically much smaller than the bifurcation 
value. 

3.2. Protostellar collapse 

What is the further fate of the clumps discussed 
above? In particular, one would like to test the no- 
tion that such a Jeans unstable clump is the imme- 
diate progenitor of a single star which forms in its 
center. That will evidently only be correct if the 
clump does not undergo further subfragmentation 
upon collapsing to higher densities. It has long 
been suspected that such subfragmcntation could 
occur at densities in excess of ~ 10 8 cm -3 , at which 
point three-body reactions become very efficient 
in converting the atomic gas (with only a trace 
amount of H2 molecules from the H~ channel) into 
fully molecular form: 3H — > H2 + H (Palla et al. 
1983). As the H2 coolant is now suddenly more 
abundant by a factor of ~ 10 3 , the corresponding 
boost in cooling could trigger a thermal instabil- 
ity, thus breaking up the clump into smaller pieces 



7 




B 1C 12 14 16 18 a 10 1Z 14 16 IF 

Lo£(r(cm)) leg[r[cm]> 



Fig. 4. — Collapse of a primordial protostar. 
Shown is the evolutionary sequence as obtained 
in a ID calculation, with time progressing from 
curves labeled 1 to 7. (a) Number density vs. ra- 
dial distance, (b) Gas temperature vs. distance, 
(c) Radial velocity vs. distance, (d) Hydrogen 
molecule abundance vs. distance. The curves la- 
beled 6 correspond to the situation briefly after 
the formation of a central hydrostatic core. (From 
Omukai & Nishi 1998.) 

(Silk 1983). Both ABN and BCL have included 
the three-body reactions in their chemical reac- 
tion networks, and have followed their simulations 
to higher densities to test whether subfragmenta- 
tion does occur or not. Both groups report that no 
further subfragmcntation is seen. With hindsight, 
that may not be too surprising. The reason be- 
ing that any small density fluctuations which are 
present earlier on, and which could serve as seeds 
for later fragmentation, will have been erased by 
pressure forces during the slow, quasi-hydrostatic 
'loitering' phase at n <~ n c . In addition, inefficient 
cooling may also play a role in suppressing high- 
density fragmentation. Despite the increase in the 
cooling rate throughout the fully molecular gas, 
this never leads to a significant drop in tempera- 
ture due to the countervailing effect of compres- 
sional heating. 



Extending the analogous calculation for the col- 
lapse of a present-day protostar (Larson 1969) to 
the primordial case, Omukai & Nishi (1998) have 
carried out one-dimensional hydrodynamical sim- 
ulations in spherical symmetry. They also consider 
the full set of chemical reactions, and implement 
an algorithm to solve for the radiative transfer in 
the H 2 lines, as well as in the continuum. The 
most important result is that the mass of the ini- 
tial hydrostatic core, formed in the center of the 
collapsing cloud when the density is high enough 
(n ~ 10 22 cm~ 3 ) for the gas to become optically 
thick to continuum radiation, is almost the same 
as the (second) core in present-day star formation: 
M core ~ 5x 10~ 3 M Q . In Figure 4, the radial 
profiles of density, temperature, velocity, and H2 
abundance are shown (reproduced from Omukai & 
Nishi 1998). The profiles of density and velocity 
before the time of core formation (corresponding 
to the curves with label 6) are well described by 
the Larson-Penston (LP) similarity solution. Once 
the core has formed, the self-similarity is broken. 
Similar results have been found by Ripamonti et 
al. (2002) who have in addition worked out the 
spectrum of the radiation that escapes from the 
collapsing clump (mostly in the IR, both as con- 
tinuum and line photons). This type of approx- 
imately self-similar behavior seems to be a very 
generic result of collapse with a simple equation 
of state, even when rotation and magnetic fields 
are included (see Larson 2003). The apparent ro- 
bustness of the LP solution thus supports the re- 
sults found for the Population III case. The small 
value of the initial core does not mean that this 
will be the final mass of a Population III star. This 
instead will be determined by how efficient the ac- 
cretion process will be in incorporating the clump 
mass into the growing protostar. 

3.3. Accretion Physics 

How massive were the first stars? Star for- 
mation typically proceeds from the 'inside-out', 
through the accretion of gas onto a central hy- 
drostatic core. Whereas the initial mass of the 
hydrostatic core is very similar for primordial and 
present-day star formation (see above), the accre- 
tion process - ultimately responsible for setting 
the final stellar mass - is expected to be rather 
different. On dimensional grounds, the accretion 
rate is simply related to the sound speed cubed 
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Fig. 5. — Evolution of accreting metal- free proto- 
star. Shown is the radius-mass relation for differ- 
ent values of the accretion rate (increasing from 
bottom to top). Accretion is effectively shut off 
for the cases with M > M cr it because of the dra- 
matic increase in radius. (From Omukai & Palla 
2003.) 

over Newton's constant (or equivalently given by 
the ratio of the Jeans mass and the free-fall time): 
Ms_ cc ~ c 3 /G oc T 3 / 2 . A simple comparison of 
the temperatures in present-day star forming re- 
gions (T ~ 10 K) with those in primordial ones 
(T ~ 200 — 300 K) already indicates a difference 
in the accretion rate of more than two orders of 
magnitude. 

Recently, Omukai & Palla (2001, 2003) have in- 
vestigated the accretion problem in considerable 
detail, going beyond the simple dimensional argu- 
ment given above. Their computational technique 
approximates the time evolution by considering 
a sequence of steady-state accretion flows onto a 
growing hydrostatic core. Somewhat counterintu- 
itively, these authors identify a critical accretion 
rate, Merit ~4x 10~ 3 M Q yr _1 , such that for ac- 
cretion rates higher than this, the protostar cannot 
grow to masses much in excess of a few 10M Q . For 
smaller rates, however, the accretion is predicted 
to proceed all the way up to ~ 600M©, i.e., of 
order the host clump. 

The physical basis for the critical accretion 
rate is that for ongoing accretion onto the core, 
the luminosity must not exceed the Eddington 
value, Ledd- In the early stages of accretion, 
before the onset of hydrogen burning, the lumi- 



nosity is approximately given by L to t ~ Luce — 
GM*M acc /.R*. By demanding L acc ~ L E dd it fol- 
lows 

Merit -^f^~5xlO- 3 M yr-\ (2) 

where i?„ ~ 5i? Q , a typical value for a Popula- 
tion III main-sequence star (e.g., Bromm et al. 
2001b). In Figure 5 (from Omukai & Palla 2003), 
the mass-radius relation is shown for various val- 
ues of the accretion rate. As can be seen, the 
dramatic swelling in radius effectively shuts off ac- 
cretion at M* < 100M©, when the accretion rate 
exceeds M cr ; t . 

Realistic accretion flows are expected to have 
a time-dependent rate, and the outcome will thus 
depend on whether the accretion rate will decline 
rapidly enough to avoid exceeding the Edding- 
ton luminosity at some stage during the evolution. 
The biggest caveat concerning the Omukai & Palla 
results seems to be the issue of geometry. A three- 
dimensional accretion flow of gas with some resid- 
ual degree of angular momentum will deviate from 
spherical symmetry, and instead form a disk. It 
is then conceivable that most of the photons can 
escape along the axes whereas mass can flow in 
unimpeded through the accretion disk (see Tan & 
McKee 2003). 

As described above, the accretion process may 
be able to incorporate a large part of the parent 
clump into the central star. Can a Population III 
star ever reach this asymptotic mass limit? The 
answer to this question is not yet known with any 
certainty, and it depends on whether the accretion 
from the dust-free envelope is eventually termi- 
nated by feedback from the star (e.g., Omukai & 
Palla 2001, 2003; Omukai & Inutsuka 2002; Ripa- 
monti ct al. 2002; Tan & McKee 2003). The stan- 
dard mechanism by which accretion may be ter- 
minated in metal-rich gas, namely radiation pres- 
sure on dust grains (Wolfire & Cassinelli 1987), 
is evidently not effective for gas with a primordial 
composition. Recently, it has been speculated that 
accretion could instead be turned off through the 
formation of an H II region (Omukai & Inutsuka 
2002), or through the radiation pressure exerted 
by trapped Lya photons (Tan & McKee 2003). 
The termination of the accretion process defines 
the current unsolved frontier in studies of Popula- 
tion III star formation. 
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3.4. The Second Generation of Stars: 
Critical Metallicity 

How and when did the transition take place 
from the early formation of massive stars to that of 
low-mass stars at later times? The second genera- 
tion of stars formed under conditions that were 
much more complicated again due to the feed- 
back from the first stars on the IGM, both due 
to the production of photons and of heavy el- 
ements. In contrast to the formation mode of 
massive stars (Population III) at high rcdshifts, 
fragmentation is observed to favor stars below a 
solar mass (Population I and II) in the present- 
day universe. The transition between these fun- 
damental modes is expected to be mainly driven 
by the progressive enrichment of the cosmic gas 
with heavy elements, which enables the gas to cool 
to lower temperatures. The concept of a 'critical 
metallicity', Z CT [ t , has been used to characterize 
the transition between Population III and Popu- 
lation II formation modes, where Z denotes the 
mass fraction contributed by all heavy elements 
(Omukai 2000; Bromm ct al. 2001a; Schneider 
ct al. 2002, 2003a; Mackey et al. 2003). These 
studies have constrained this important param- 
eter to only within a few orders of magnitude, 
-^crit ~ 10~ 6 — 10~ 3 Zq, under the implicit assump- 
tion of solar relative abundances of metals. This 
assumption is likely to be violated by the metal 
yields of the first SNe at high-redshifts, for which 
strong deviations from solar abundance ratios are 
predicted (e.g., Oh et al. 2001; Hcger & Woosley 
2002; Qian, Sargent, & Wasserburg 2002; Qian & 
Wasserburg 2002; Umeda & Nomoto 2002, 2003). 

Recently, Bromm & Loeb (2003b) have shown 
that the transition between the above star forma- 
tion modes is driven primarily by fine-structure 
line cooling of singly-ionized carbon or neu- 
tral atomic oxygen. Earlier estimates of Z CT - lt 
which did not explicitly distinguish between dif- 
ferent coolants are refined by introducing sep- 
arate critical abundances for carbon and oxy- 
gen, [C/H] cr ;t and [0/H] cr ; t , respectively, where 
[A/H]= log 10 (iVA/^H) - \og w (N K /N n ) Q . Since 
C and O are also the most important coolants 
throughout most of the cool atomic ISM in 
present-day galaxies, it is not implausible that 
these species might be responsible for the global 
shift in the star formation mode. Under the tem- 
perature and density conditions that characterize 



Population III star formation, the fine-structure 
lines of O I and C II dominate over all other 
metal transitions (see Hollenbach & McKee 1989). 
Cooling due to molecules becomes important only 
at lower temperatures, and cooling due to dust 
grains only at higher densities (e.g., Omukai 2000; 
Schneider ct al. 2003a). The presence of dust 
is likely to modify the equation of state at these 
high densities in important ways (Schneider et al. 
2003a), and it will be interesting to explore its 
role in future collapse calculations. The physi- 
cal nature of the dust, however, that is produced 
by the first SNe is currently still quite uncer- 
tain (e.g., Loeb & Haiman 1997; Todini & Fer- 
rara 2001; Nozawa ct al. 2003; Schneider, Fer- 
rara, & Salvaterra 2003b). Numerically, the crit- 
ical C and O abundances are estimated to be: 
[C/H] crit ^ -3.5 ±0.1 and [0/H] crit ~ -3.1 ±0.2. 

Even if sufficient C or O atoms are present to 
further cool the gas, there will be a minimum at- 
tainable temperature that is set by the interac- 
tion of the atoms with the thermal CMB: Tcmb = 
2.7 K(l + z) (e.g., Larson 1998; Clarke & Bromm 
2003). At z ~ 15, this results in a characteris- 
tic stellar mass of M* ~ 20M Q (n / /10 4 cm" 3 )- 1 / 2 , 
where n/ > 10 4 cm -3 is the density at which 
opacity prevents further fragmentation (e.g., Rees 
1976). It is possible that the transition from 
the high-mass to the low-mass star formation 
mode was modulated by the CMB temperature 
and was therefore gradual, involving intermediate- 
mass ('Population II. 5') stars at intermediate red- 
shifts (Mackey et al. 2003). This transitional pop- 
ulation could give rise to the faint SNe that have 
been proposed to explain the observed abundance 
patterns in metal-poor stars (Umeda & Nomoto 
2002, 2003). When and how uniformly the tran- 
sition in the cosmic star formation mode did take 
place was governed by the detailed enrichment his- 
tory of the IGM. This in turn was determined by 
the hydrodynamical transport and mixing of met- 
als from the first SN explosions (e.g., Mori et al. 
2002; Bromm et al. 2003; Scannapieco et al. 2003; 
Wada & Venkatesan 2003). The transport and the 
mixing of the first heavy elements into the pristine 
IGM is currently another subject at the frontier of 
astrophysical cosmology. 
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3.5. Stellar Archaeology: Relics from the 
End of the Dark Ages 

It has long been realized that the most metal- 
poor stars found in our cosmic neighborhood 
would encode the signature from the first stars 
within their elemental abundance pattern (e.g., 
Bond 1981; Beers, Preston, & Shectman 1992). 
For many decades, however, the observational 
search has failed to discover a truly first-generation 
star with zero metallicity. Indeed, there seemed to 
have been an observational lower limit of [Fc/H] 
<~ —4 (e.g., Carr 1987). In view of the recent 
theoretical prediction that most Population III 
stars were very massive, with associated lifetimes 
of ~ 10 6 yr, the failure to find any 'living' Popu- 
lation III star in the Galaxy is not surprising, as 
they would all have died a long time ago (e.g., Her- 
nandez & Ferrara 2001). Furthermore, theory has 
predicted that star formation out of extremely 
low-metallicity gas, with Z < Z alt ~ 10~ 3 - 5 Z Q 
(see § 3.4), would be essentially equivalent to that 
out of truly primordial gas. Again, this theoreti- 
cal prediction was in accordance with the apparent 
observed lower-metallicity cutoff. 

Recently, however, this simple picture has been 
challenged by the discovery of the star HE0107- 
5240 with a mass of O.8M and an iron abundance 
of [Fc/H] = -5.3 (Christlieb et al. 2002). This 
finding indicates that at least some low mass stars 
could have formed out of extremely low-metallicity 
gas. Does the existence of this star invalidate the 
theory of a metallicity threshold for enabling low- 
mass star formation? As pointed out by Umeda 
& Nomoto (2003), a possible explanation could lie 
in the unusually high abundances of carbon and 
oxygen in HE0107-5240. 

In Figure 6, the theoretical C and O thresh- 
olds derived by Bromm & Loeb (2003b) are com- 
pared to the observed abundances in metal-poor 
dwarf (Akcrman et al. 2003) and giant (Cayrel 
et al. 2003) stars in the halo of our Galaxy. As 
can be seen, all data points lie above the critical 
O abundance but a few cases lie below the criti- 
cal C threshold. All of these low mass stars are 
consistent with the model since the correspond- 
ing O abundances lie above the predicted thresh- 
old. The sub-critical [C/H] abundances could have 
either originated in the progenitor cloud or from 
the mixing of CNO-processed material (with car- 
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Fig. 6. — Observed abundances in low-metallicity 
Galactic halo stars. For both carbon (upper panel) 
and oxygen (lower panel), filled circles correspond 
to samples of dwarf and subgiant stars (from Ak- 
crman et al. 2003), and open squares to a sam- 
ple of giant stars (from Cayrel et al. 2003). The 
dashed lines indicate the predicted critical carbon 
and oxygen abundances (see § 3.4). Highlighted 
(cross) is the location of the extremely iron-poor 
giant star HE0107-5240. (Adapted from Bromm 
& Loeb 2003b.) 



bon converted into nitrogen) into the stellar atmo- 
sphere during the red giant phase. Note that the 
extremely iron-poor star HE0 107-5240 has C and 
O abundances that both lie above the respective 
critical levels. The formation of this low mass star 
(~ 0.8M Q ) is therefore consistent with the theo- 
retical framework considered by Bromm & Loeb 
(2003b). 

The lessons from stellar archaeology on the na- 
ture of the first stars are likely to increase in im- 
portance, since greatly improved, large surveys of 
metal-poor Galactic halo stars are under way, or 
are currently being planned. 
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4. THE FIRST QUASARS 

Quasars are believed to be powered by the heat 
generated during the accretion of gas onto SMBHs 
(e.g., Rees 1984). The existence of SMBHs with 
inferred masses of > 10 9 M©, less than a billion 
years after the big bang, as implied by the recent 
discovery of quasars at redshifs z > 6 (e.g., Fan et 
al. 2003), provides important constraints on the 
SMBH formation scenario (Haiman & Loeb 2001). 
Can the seeds of SMBHs form through the direct 
collapse of primordial gas clouds at high redshifts? 
Previous work (Loeb & Rasio 1994) has shown 
that without a pre-existing central point mass, this 
is rendered difficult by the negative feedback re- 
sulting from star formation in the collapsing cloud. 
The input of kinetic energy due to supernova ex- 
plosions prevents the gas from assembling in the 
center of the dark matter potential well, thus pre- 
cluding the direct formation of an SMBH. If, how- 
ever, star formation was suppressed in a cloud that 
could still undergo overall collapse, such an ad- 
verse feedback would not occur. 

The SMBH formation problem has recently 
been revisited with SPH simulations of isolated 
2(T— peaks with total masses of 10 8 M Q that col- 
lapse at z v i r ~ 10 (Bromm & Loeb 2003a). The 
virial temperature of these dwarf galaxies exceeds 
~ 10 4 K, so as to allow collapse of their gas through 
cooling by atomic hydrogen transitions (Oh & 
Haiman 2002). Since structure formation proceeds 
in a bottom-up fashion, such a system would en- 
compass lower-mass halos that would have col- 
lapsed earlier on. These sub-systems have virial 
temperatures below 10 4 K, and consequently rely 
on the presence of H2 for their cooling. Molecular 
hydrogen, however, is fragile and readily destroyed 
by photons in the Lyman- Werner bands with en- 
ergies (11.2 — 13.6 eV) just below the Lyman limit 
(Haiman, Rees, & Loeb 1997). These photons are 
able to penetrate a predominantly neutral IGM. 

At first, the limiting case is considered where H2 
destruction is complete. Depending on the initial 
spin, either one (for zero initial spin) or two com- 
pact objects form with masses in excess of 10 6 M Q 
and radii < 1 pc (see Fig. 7). In the case of nonzero 
spin a binary system of clumps has formed with 
a separation of ~ 10 pc. Such a system of two 
compact objects is expected to efficiently radiate 
gravitational waves that could be detectable with 
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Fig. 7. — Central gas density in a dwarf galaxy 
with virial temperature just above the atomic cool- 
ing threshold, but with no H2 molecules present. 
Shown is the projection in the x — y plane at 
z ~ 10. The box size is 200 pc on a side. In this 
case, the initial spin is A = 0.05. Two compact 
objects have formed with masses of 2.2 x 10 6 M Q 
and 3.1 x 10 6 M©, respectively, and radii < 1 pc. 
(Adapted from Bromm & Loeb 2003a). 

the planned Laser Interferometer Space Antenna 2 
(LISA; Wyithe & Loeb 2003c). 

What is the further fate of the central ob- 
ject? Once the gas has collapsed to densities 
above ~ 10 17 cm -3 and radii < 10 16 cm, Thom- 
son scattering traps the photons, and the cooling 
time consequently becomes much larger than both 
the free-fall and viscous timescales (see Bromm 
& Loeb 2003a for details). The gas is therefore 
likely to settle into a radiation-pressure supported 
configuration resembling a rotating supermassive 
star. Recent fully-relativistic calculations of the 
evolution of such stars predict that they would in- 
evitably collapse to a massive black hole (Baum- 
garte & Shapiro 1999; Shibata & Shapiro 2002). 
Under a wide range of initial conditions, a sub- 
stantial fraction (~ 90%) of the mass of the su- 
permassive star is expected to end up in the black 
hole. 



2 See http://lisa.jpl.nasa.gov/ 
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Is such a complete destruction of H2 possible? 
When including an external background of soft 
UV radiation in the simulation, Bromm & Loeb 
(2003a) find that a flux level comparable to what 
is expected close to the end of the reionization 
epoch is sufficient to suppress H 2 molecule forma- 
tion. This is the case even when the effect of self- 
shielding is taken into account (Draine & Bertoldi 
1996). The effective suppression of H 2 formation 
crucially depends on the presence of a stellar-like 
radiation background. It is therefore likely that 
stars existed before the first quasars could have 
formed. 

5. CONCLUSIONS 

The first stars and quasars, their formation and 
their cosmological implications, are a fascinating 
subject at the frontier of modern cosmology. Up 
to now, most of our efforts have been theoretical, 
but we are approaching the point where observa- 
tions can test our ideas. This is very significant, as 
we are bound to learn important lessons about the 
physical state of the universe at the end of its dark 
age. These empirical probes are provided by both 
high-redshift observations, and complementary to 
this, by the 'near-field cosmology' (Freeman & 
Bland-Hawthorn 2002) of scrutinizing the chem- 
ical abundance patterns in extremely metal-poor 
stars in our local cosmic neighborhood. One can 
only wonder what exciting discoveries are awaiting 
us, and it is a great privilege to be a part of this 
grand endeavor. 
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